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^^ ' Abstract: The short distance behavior of baryon-baryon potentials defined through 

Nambu-Bethe-Salpeter wave functions is investigated using the operator product expan- 
sion. In a previous analysis of the nucleon-nucleon case, corresponding to the SU(3) chan- 

k> , nels 27s and 10a, we argued that the potentials have a repulsive core. A new feature occurs 

H ' for the case of baryons made up of three flavors: manifestly asymptotically attractive po- 

tentials appear in the singlet and octet channels. Attraction in the singlet channel was 
first indicated by quark model considerations, and recently been found in numerical lattice 
simulations. The latter have however not yet revealed asymptotic attraction in the octet 
channels; we give a speculative explanation for this apparent discrepancy. 
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1. Introduction 

In a recent paper Q a proposal has been made to study nucleon-nucleon (NN) potentials 
from a first principle QCD approach. In this field theoretic framework, potentials are 
obtained through the Schrodinger operator applied to Nambu-Bethe-Salpeter (NBS) wave 
functions: 

v(j) = E+iZW)^ (1.1) 

2id (Pe{x) 
Here ^ is the two nucleon reduced mass and the NBS wave function is defined by 

ipEix) = {0\N{x/2, 0)N{-x/2, 0) |2N, E) , (1.2) 

where |2N, ii^) is a QCD eigenstate with energy E (suppressing here other quantum num- 
bers), and A^ is a nucleon interpolating operator made of 3 quarks. Such wave func- 
tions have been measured through numerical simulations of the lattice regularized theory 
m, §, m, ^. Although many conceptual questions remain to be resolved (the potential 



defined by (LI) ^^ ^^^ example energy dependent), the corresponding potentials indeed 
qualitatively resemble phenomenological NN potentials which are widely used in nuclear 
physics. The force at medium to long distance (r > 2 fm) is shown to be attractive. This 
feature has long well been understood in terms of pion and other heavier meson exchanges. 
At short distance, a characteristic repulsive core is reproduced by the lattice QCD sim- 
ulation ffl. In P] we performed an operator product expansion (OPE) analysis of NN 
NBS wave functions in QCD to theoretically better understand the repulsive core of the 
NN potential (at least that of the measured NBS potential). Thanks to the property of 
asymptotic freedom of QCD the form of leading short distance behavior of the coefficent 
functions can be computed using perturbation theory. To set the stage we first give a very 
short summary of the results of fsl . 



The behavior of the wave functions (Pe{x) at short distances (r = \x\ — )• 0) is encoded 
in the operator product expansion (OPE) of the two nucleon operators: 

N{x/2,0)N{-x/2,0)^Y.^k{x)Ok{0,0), (1.3) 

k 

where {O^} is a set of local color singlet 6-quark operators with two-nucleon quantum 
numbers. Asymptotically the x-dependence and energy dependence of the wave function 
is factorized into 

ifEix) ^Y.^k{x){0\Oki6,0)\2N,E) . (1.4) 

A; 

Standard renormalization group (RG) analysis gives Q the leading short distance behavior 
of the OPE coefficient function as 

Dk{x) ^ {in'^y dk , (1.5) 

where Uk is related to the 1-loop coefficient of the anomalous dimension of the operator Ok, 
dk is the tree-level contribution^ of -Dfc(O) and finally ro is some typical non-perturbative 
QCD scale (which we take here to be 1 fm). Clearly the operator with largest RG power v^ 
dominates the wave function ( [1.2| ) at short distances. We will denote the leading (largest) 
power by ui and the subleading one (second largest) by 1^2 ■ 

If neither vi nor di vanishes, this leads to the leading asymptotics of the s-wave po- 
tential of the form 

V{r)^- J\^^ , (1.6) 



r 



(In^) 



which is attractive for z^i > and repulsive for ui < 0. Note that the above asymptotic 
form is energy independent. 

If vi = 0, the situation is more complicated. In this case the relative sign of the ratio 
R between the leading contribution and the subleading contribution, corresponding to an 
operator O2 with z^2 < 0, is important. If R is positive, the potential is repulsive, while it 
is attractive for negative R. The leading asymptotics of the potential is energy dependent 
in this case, since it is proportional to R which depends on the energy and it may even 
change sign as function of E. 

In the nucleon potential problem the above degenerate case {ui = 0) occurs. In |5| we 
have argued that the relative coefficient R is positive and the short distance limit of the 
nucleon potential is repulsive in both possible isospin-spin s-wave channels {IJ = 01 and 
/J = 10). 

To actually compute the renormalization group power u^ of an operator, we first have 
to compute the 1-loop coefficient of its anomalous dimension defined by [g] : 

7fc(5)=7i'^5' + 0(/). (1.7) 



^In principle, the leading term might come from an operator that has vanishing tree level coefficient 
dfc and enters the two-nucleon OPE at higher order in perturbation theory, provided that its anomalous 
dimension is large enough to compensate. In |5| we calculated the anomalous dimensions of all relevant 
6-quark operators and showed that no such operators exist. 



Here we assume that the corresponding operator O^ is renormahzed multiphcatively. It is 
always possible to choose an operator basis where this is the case, at least up to 1-loop 
level. Having computed the coefficient 7^ the RG power v^ is given by 

^ 7l"-27!J' 48.^71" -24 

where j3q is the 1-loop beta function coefficient of the Nf flavor theory and 7]^ = 12/487r^ 
is the 1-loop coefficient of the anomalous dimension of the nucleon operator. For simplicity 
we count the anomalous dimension coefficients in units of l/487r^ and for example we shall 
call an operator O^ of anomalous dimension 36 if 7^ = j^ = -^^ . 

In this paper we extend the analysis of ref. ||5[ to the case of three flavors. Instead of 
nucleons, we shall consider the baryon octet and study the short distance behavior of the 
NBS wave function 

^%iS) = (0|i?i(x/2,0)i?^(-i;/2,0)|2B,i?), (1.9) 

where i,j are SU(3) flavor octet indices, a, /3 spin indices, S^ baryon field operators and 
|2B, £^) 2-baryon states with energy E (and suppressed other quantum numbers). Since the 
operator anomalous dimensions are mass-independent in perturbation theory, all our results 
are manifestly SU(3) symmetric. The tensor product of two octets can be decomposed as 

8^8 = (1 080 27), + (80 100 IO)a, (1.10) 

where the first three representations are symmetric, the last three antisymmetric in the 
two baryon indices. In terms of octet fields, for example, we have 

BB^^''\l = 0,13 = 0,y = 0) = y^AA - y^SS + y^iVH, (1.11) 

i?S(«^) (0,0,0) = -yiAA-y|sS + yiiVH, (1.12) 

55^(0,0,0) = -W^AA + A/lsS + A/^iVH, (1.13) 



with 



and 



SS = y ^S+S- - y ^S^SO + y ^S-S+ (1.14) 

iVH = y|p- + yis-p - yin-° - yi-On (1.15) 



Bi?(i°)(0,0,2) = y-pn-y-np, (1.16) 

i?S(io) (3/2, 3/2,1) = y|pS+ - y|s+p, (1.17) 

BS^S'^) (0,0,0) = yipH- - y|H-p - Tins" + yisOn . (1.18) 



Here / and I3 are the total isospin and its third component while Y is the total hyper charge 
for two baryons. The spin quantum numbers are understood as in ( |2.5| ). There are six 
s-wave potentials corresponding to these six channels, since, due to the total antisymmetry 
of the two-baryon product, the flavor representation and its symmetry properties uniquely 
determine the spin representation for s-waves: with the antisymmetric representations the 
symmetric J = 1 spin representation is coupled and J = is paired with the symmetric 
flavor representations. In this scheme the nucleon potentials are in the 27 and 10 channels, 
for isospin-spin IJ = 10 and IJ = 01, respectively. 

2. 3-flavor results 

We have repeated the analysis of Q for the case of three quark flavors. In the three-flavor 
case some channels may become attractive at short distance since the Pauli exclusion 
principle is less significant than in the 2-flavor case. In this paper we therefore investigate 
which representations of the flavor SU(3) become attractive at short distance using the 
OPE. The method and the formulae necessary to compute the 1-loop anomalous dimension 
of any local gauge invariant 6-quark operator are given in detail in ref. [^ and will not be 
reproduced here. 

The local gauge invariant 6-quark operators we need are local products of six quark field 
operators qa with flavor index / € {u, d, s}, color index a and Dirac index a E {1, 2, 3, 4} 
and can be written as local products 

rpfah rpf'g'h' /r, i N 

of two color singlet triple-quark objects of the form 

Tit, = e^'" q^QlV," = fa9^h, (2.2) 

where f,g,h,f',g',h' are flavor indices and a, I3,^,a' , P' ,^' Dirac indices. The Dirac matri- 
ces are chosen such that a = 1,2 correspond to positive chirality and a = 3,4 to negative 
chirality. As shown in |5|, the 1-loop renormalization matrix mixes operators with one 
pair of flavor and/or Dirac indices interchanged, provided that the Dirac indices belong to 
the same chirality. From this it follows that operators with different flavor or Dirac index 
structure do not mix. It also follows that (for any flavor structure) we always have oper- 
ators with vanishing RG power since for the Dirac structure (the square brackets indicate 
antisymmetrization) 

'■1[12] ^3 [34] 

index exchange between quark fields belonging to different triplets is not allowed by the 
above mentioned chirality rule and therefore the operator has anomalous dimension 24 
(coming from exchange within the triplets), which is exactly compensated by the contri- 
bution of the two baryons as in ( |l.8|) . 

It is sufficient to consider 6-quark operators with flavor structure uuddss since all six 
channels of ( 1.10| ) are present in this case. There are many possible Dirac structures but 



rpfgh rpf'g'h' /r, o\ 



there are only a few among these that contain "attractive" operators (i.e. operators with 



anomalous dimension > 24 in our units). We have considered all Dirac structures, calcu- 
lated the 1-loop anomalous dimension matrix and performed the change of basis necessary 
to make operators renormalize multiplicatively (diagonally) at 1-loop level. We refrain 
from giving a complete list here (which would be quite long: with Dirac structure 112334 
for example there are 123 independent uuddss flavor 6-quark operators), and give here the 
complete list of attractive operators only. These are 

• Dirac structure: 112334, anomalous dimension: 42 (SU(3) singlet) 

• Dirac structure: 111222, anomalous dimension: 36 (SU(3) singlet) 

• Dirac structure: 112234, anomalous dimension: 36 (SU(3) singlet and octet) 

• Dirac structure: 111234, anomalous dimension: 32 (SU(3) singlet and octet) 

• Dirac structure: 112234, anomalous dimension: 32 (SU(3) singlet and octet) 

We here listed operators with chiral (1,2) indices not less in number than antichiral (3,4) 
ones. There are also parity reflected operators with the same anomalous dimension and 
flavor quantum numbers. 

We still have to study the question whether these operators are already present at tree 
level in the OPE since otherwise their effective anomalous dimension is reduced by a large 
number (54 in our units for every extra order in perturbation theory with three flavors) 
and they become repulsive. 

The baryon field operators corresponding to the flavor structure fgh are built from 
the building blocks 

^aL = fa{gih2 " 92^) , b{f^ = faig^h^ - Qih^i) , (2.4) 

for a = 1, 2, 3, 4 and the physical baryon field is 

where a = 1,2 and 1 = 3, 2 = 4. We sometimes use octet indices {p = uud, E~ = ssd, 
etc.) instead of the flavor notation fgh. 

It turns out that (with the exception of the anomalous dimension 32 flavor singlet 
operators) all the attractive operators are indeed present at tree level in the OPE for two 
building blocks of the form ( |2.4D . We find 

bl.L% + bl.nbl, - 5^^-7t''^'li^lf + •.., (2.6) 

+ d*J'^ /3[36;8] e^^ ^[36]fc ^ _ _ _ ^ (2.7) 

b^L^^^ ^ 5*^ap''^^]e„;3M[36] + ..., (2.8) 

hl^hl ^ <5*^af^^]e„;3if'^]+d^^''^a['''^'U„;,L[36]^ + ..., (2.9) 



where . . . stand for non-attractive operators. We also have the parity reflected counterparts 
of the above relations. We will denote the parity reflected operators with a . Here the 
coefficients 7, /3, a are nonvanishing and the K, L, M are local 6-quark operators with their 
anomalous dimension and SU(3) quantum numbers indicated. 
Finally we can compute the baryon-baryon tree level OPE: 

Bl^Bj « 6^^ e^p {a,Af^ + a,Af^ + a,Af^ } + d^^^ e^^ bB^'^^'' + P^' cC^S^' + . . . (2.10) 

Here 

a, = af'-^'\ 436] ^^[36] ^^^[36]^ 

as = ^^''-''\ 4^1=i,Jf-<l, (2.11) 

b = a[36;8] _^ ^[36;8] ^ g[36]fc ^ ^[36]fc _^ ^[36]fc^ 

_ /o[32;8] ^[32]fc _ . [32]fe ? [32]fc 

'' ~ P ' '-a/3 - -'^a/J + ^a/3 • 

We see that there are no attractive operators in the 27^, 10a and 10a channels. This is 
consistent with what we found in the nucleon case; these belong to 27^ and 10a • There are 
three attractive channels: 1^, 8^, 8a. Among these the NBS wave function for Ig becomes 
most singular at short distance since its anomalous dimension 42 is largest. In terms of the 
corresponding potential, however, all three have the same r dependence at short distance. 
They may not be equally strong: the coefficients ai etc. can be calculated^ but these 
coefficients are to be multiplied by the matrix elements of the corresponding operators, 
which we do not know. 

We note that there is no local 6-quark operator in the 8^ channel with nonvanishing 
non-relativistic limit, and therefore the corresponding matrix elements also vanish in the 
same limit. We considered a class of 6-quark operators in the non-relativistic limit, which 
have uuddss flavor and Sz = with Sz being the z component of the total spin. Since qa = 
qSi in the non-relativistic limit, this class of operators has 111222 Dirac structure. We found 
by an explicit calculation that all attractive 6-quark operators in the 8s channel vanish in 
the non-relativistic limit. Next we calculated the eigenvalues T of the SU(3) Casimir 
operator (and their eigenvectors) for this class of operators. The octet representation 
corresponds to T = 3. We have found that while the operator piE^ — H^j?2 contains T = 3 
components, the operator piH^ -|- E^p2 does not contain T = 3 components at all. (The 
subscripts 1,2 represent spin indices.) This shows that there is no 8^ 6-quark operator in 
the non-relativistic limit. 

One can also see the above property group theoretically [^J. The 6-quark operator must 
be totally antisymmetric due to the Fermi-Dirac statistics. In terms of the permutation 
group, the totally antisymmetric representation corresponds to [1^], whose Young diagram 
has 6 boxes in a column. Each quark operator in the non-relativistic limit has color (c), 
spin(s), flavor (/) in addition to the coordinate x which is common for all 6 quarks in the 



^e.g. at tree level a[36.8]/^[36,8] ^ .g/g ^nd of •^V/?"'''' = 1/3 



local operator. We now consider symmetries of each index. The symmetry for color SU(3) 
must be [2^]c due to gauge invariance, while that for the coordinate is totally symmetric, 
therefore [6]^;. For 3 flavor, the non-relativistic quark belongs to the fundamental repre- 
sentation of SU(6)s/ D SU(2)s® SU(3)j. The symmetry of 6 quarks in SU(6)sj, denoted 
as [F]sf, must satisfy 



[2% ^ [F]sf [6]x D [11 ^ [21c (8) [Flf D [l^, (2.12) 

in order to make the 6-quark operator totally antisymmetric. In SU(6)sj, [3] contains octet 
and decuplet, and therefore two octet baryon state belongs to [3] (g) [3] = [6] © [4, 2] © [5, 1] © 
[3^], where the first(last) two are (anti)symmetric under exchange of two [3]. Since two 
baryons at the same space-time point must be antisymmetric under exchange of the two, 
[5, 1] and [3^] are allowed. The table C on p. 116 of Ref.|0 tells us 

[2% (S) [5, 1],/ = [2^, l2] e [3, 2, 1], [2^], [3^],^ = [l6] © [22, l^] © [3^] © [4, l2].(2.13) 



These show that the [5, 1]^^ representation does not satisfy the condition ( |2.12| ) while ['i'^]sf 
does. Therefore only [32]^^, but not [5,1]^/, must appear in the local (gauge-invariant) 6- 
quark operator. In terms of SU(6)s/, the 6 quarks in the 8^ channel, which has [3^]^ in 
spin and [3, 2, l]f in flavor, becomes (see also the same table C) 

8s = [3% ® [3, 2, 1]/ = [5, 1] © [4, 2] © [4, 1^] © 2[3, 2, 1] © 3[3, 1^] © [2^, 1^] e [2, 1^] 

(2.14) 



while Is and 8a are expressed as 

If 



Is = [3% © [2^]; = [3^] © [4, l2] © [22, l2] e [l6] (2.15) 



8a = [4, 2], © [3, 2, l]f = [5, 1] © 2[4, 2] © 2[4, 1^] © [3^] © 3[3, 2, 1] 

©[2^] © 2[3, 1^] © 2[22, l2] e [2, 1^]. (2.16) 

Since 8^ contains [5, 1] only but no [3^], the 6 quark operator in this channel does not exist 
in the 811(6)^/ representation (i.e. the non-relativistic limit). 

3. Comparison to quark model calculations and lattice Monte Carlo 

The short distance part of the baryon potential has also been calculated using a valence 
quark model with gluon exchange [^ and lattice Monte Carlo measurements [§]. 

In the quark model above the only attractive channel is the SU(3) singlet one (which 
may correspond to a strange dibaryon bound state: the H particle). It is not possible to 
quantitatively compare the results with ours since the way the baryon potential is defined 
in the two formalisms is different and in the valence quark model SU(3) breaking effects are 
also taken into account (and they are, in some cases, comparable to the SU(3) symmetric 
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Figure 1: Potential corresponding to the model wave function with i/i — 2/9, 1/2 = —0.10 and R = 2. 



contributions). In addition, a particle basis (A^, H, S, A) is used in the valence quark model 
and the various s-wave channels are labeled as in this example: 



iVS (/,J) = (i,0). 



(3.1) 



In the quark model this channel is for example strongly repulsive but this result cannot be 
directly compared to ours since in our language this is a linear combination of the channels 
27s and 85. On the other hand, for the case 



iVS (/,J) = (-,0) 



(3.2) 



the quark model gives (a somewhat weaker) repulsion, consistently with our result, since 
this channel is pure 27^. 

The MC measurements H are performed using full QCD with three massive quarks 
in the flavor SU(3) limit. The results (for s-wave potentials) are analysed in terms of the 
six channels ( LlOf ) and attraction is found only in the singlet channel. The other five are 
found to be repulsive, with four of them (including the 27 and 10 already occurring in 
the nucleon case) being rather similar. The potential corresponding to the 8^ channel is 
different from all the other ones in that it does not have any attractive pocket and its 
repulsive core is stronger and broader than for the other cases. 

Our results do not necessarily contradict to what was found in the quark model and 
measured on the lattice since the distance range in which our asymptotic results are valid 
may not overlap with the range where the other two methods work. The latter is the 
range between 0.1 fm and 1.5 fm, whereas the true asymptotic behavior predicted by 
asymptotically free perturbation theory may be valid in a much shorter distance range 
only and in the range where comparison is possible subleading contributions in the OPE 
may be non-negligible. To illustrate this situation we consider a (renormalization group 
motivated) simpified model where the baryon wave function is the sum of a leading term 
with RG power z^i and a subleading term with relative coefficient R and RG power 1^2 • 



More precisely, we take 



Y + KlnA = ln— , K=— 2=— . (3.4) 



as our simplified baryon wave function, where A(r) is the 2-loop 3-flavor running coupling 
defined by 

A r ' 2/32 81 

Using A(r) instead of 1/ ln(ro/r) takes into account some of the higher perturbative contri- 
butions. Of course, the complete RG asymptotic form of the NBS wave function is much 
more involved than the simple formula (|3.3| ). There are perturbative corrections (power 
series in A) multiplying both the leading power of A and the subleading one and there 
are also many sub-subleading terms to be added. Nevertheless the simple model ( p. 3]) can 
qualitatively reproduce the main features of the true asymptotics of the wave function. We 
take 

1/1 = ^, 1/2 = -0.1. (3.5) 

We have chosen these numbers because we want to model the 8^ channel, which is asymp- 
totically attractive but looks strongly repulsive in the range studied by lattice Monte Carlo 
and valence quark models, ui corresponds to anomalous dimension 36, which is the largest 
(and only attractive one) occurring in this channel and 1^2 is chosen to represent the "aver- 
age" of the RG power zero (which, as we have seen, occurs in all channels) and the many 
negative sub-subleading powers. 

Clearly for r ^- eventually the leading term dominates, but since z^i and 1^2 are both 
small numbers the subleading effects are still important in the range 0.1-0.3 fm, where 
the short distance expansion (|3.3| ) can be compared to results obtained by the above two 
methods. If the relative coefficient R becomes large, the second term will dominate in this 
range and the potential becomes repulsive. In Fig. ^ the potential^ for i? = 2 is shown. It 
is less attractive than the one corresponding to the leading term alone, but is still weakly 
attractive. On the other hand for R = —2 already the second term dominates and the 
potential, which is shown in Fig. 0, is clearly repulsive in this range. 

A natural question arising is that if this mechanism explains why an asymptotically 
attractive potential can look repulsive in the intermediate range 0.1-0.3 fm, why the same 
mechanism can not modify the results obtained in the nucleon case (where repulsion was 
found in agreement with the quark model and Monte Carlo results). The answer is that in 
the nucleon case all subleading terms correspond to repulsion so if they happen to dominate 
in the intermediate range the qualitative behavior of the potential remains unchanged. 
Actually the nucleon case corresponds to the degenerate case with i^i = and the true 
asymptotics is determined by the sign of the relative coefficient R. But the above explained 
mechanism actually improves the situation in the intermediate range since if R becomes 
large (with either sign), it dominates and the potential looks repulsive. Note that the 
potential should become attractive asymptotically at short distance for negative R. The 
attractive potential can be seen in the intermediate range only if i? ~ — 1 or smaller. 



^The model potential is simply defined as U{r) = V'^ip/ip. We are using ro units for the plots. 
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Figure 2: Potential corresponding to the model wave function with ui = 2/9, ^2 ~ —0.10 and R — —2. 



In addition to the above mechanism, the matrix elements of the local 6 quark operator 
in the 8s channel vanish in the non-relativistic limit, as shown in the previous section. 
Therefore the attraction in the 8^ channel predicted by the OPE analysis can be hardly 
seen, as long as the non-relativistic approximation holds. Indeed this approximation seems 
good in the MC measurements ||8| since the wave function in the 8s channel is observed to 
be very small at the origin. 

4. Conclusion 

In this paper we extended the OPE analysis of ref. [^ from two to three flavors, in order 
to investigate the short distance behavior of general baryon-baryon potentials in flavor 
SU(3). Interestingly the OPE analysis leads to manifest attraction for the SU(3) singlet 
potential at short distance, which is consistent with the quark model prediction and the 
lattice QCD result. The manifest attraction channel does not exist in the nucleon-nucleon 
potential for 2 flavors (the channels 27s and 10^ in the 3 flavor theory). The increase of 
the number of flavors decreases the repulsion at short distance and in some cases turns 
it into attraction already for 3 flavors, since the Pauli exclusion principle among quarks 
becomes less significant for larger number of flavors. The agreement among the three 
different methods strongly corroborates that the interaction between two baryons in the 
singlet channel is indeed attractive. It would be important to check by a lattice QCD 
calculation whether this attraction in the flavor SU(3) limit leads to a bound-state (the 
H-particle). Even if the bound state is indeed present in the SU(3) limit, it is possible that 
it becomes a resonance above the AA threshold in the real world where the strange quark 
is heavier than the other two. 

On the other hand, our OPE analysis indicates that there is short distance attraction 
in both (8a and 8s) octet channels, which disagrees with present lattice QCD results. 
Although we gave a possible explanation for this discrepancy, a direct confirmation by 
lattice QCD calculations will be needed for definite conclusions. 
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